A Novel Nuclear Model for Double Beta Decay 

Franjo Krmpotic 

Instituto de Fisica, Universidade de Sao Paulo, 05315-970 Sao Paulo-SP, Brazil 
Instituto de Fisica La Plata, CONICET, 1900 La Plata, Argentina and 
Facultad de Ciencias Astronomicas y Geofisicas, 
Universidad Nacional de La Plata, 1900 La Plata, Argentina, 
(Dated: February 9, 2008) 

Abstract 

The possibility of applying the Quasiparticle Tamm-Dancoff Approximation (QTDA) to describe 
the nuclear double beta decay is explored. Several serious inconveniences found in the Quasiparticle 
Random Phase Approximation (QRPA), such as: i) the extreme sensitivity of the 2vf3f3 decay 
amplitudes M.i v on the residual interaction in the particle-particle channel, ii) the ambiguity 
in treating the intermediate states, and iii) the need for performing a second charge-conserving 
QRPA to describe the /3/3-decays to the excited final states, are not present in the QTDA. Also, 
the QTDA allows for explicit evaluation of energy distributions of the double-charge-exchange 
transition strengths and of their sum rules, and can be straightforwardly applied to single- and 
double-closed shell nuclei. As an example, the Ca— > 48 Ti decay is discussed within the 1/p-shell 
in the particle-hole limit of the QTDA. The general [(1, l)-Pade-approximant-like] behavior of 
the 2i//3/3-decay amplitude in the plain QRPA as well as within its different variations is briefly 
reviewed. 



1 



I. INTRODUCTION 



It is a great pleasure and a great honor for me to contribute to this commemorative issue 
in memory of Dubravko Tadic, who was one of my closest friends and the best coworker I 
ever had. I miss him badly, as many people do! We cooperated closely since 1966. First, we 



studied the single beta (/3)-deca 
in the double beta (/3/3)-decay 




1 ^J, and in recent years we were basically involved 
and the nonmesonic weak decay of hypernuclei 



0, 1(3, [ll|, [12I Q. These topics are nice examples of interrelation between Particle and 
Nuclear Physics. In fact, Dubravko Tadic took part in many important developments of the 
theory of weak interactions, as well as in the advancement of particle and nuclear physics as a 
whole. With the entanglement between birds and fishes in the Escher's engraving, shown in 
Figure 1, I want to symbolize the close cooperation I had with Dubravko. To tell the truth, 
I have done my first work in theoretical physics, entitled: On the Induced Terms and Partial 
Conservation of the Axial Vector Current in Beta Decay Q|, under Dubravko's guidance, and 
the line of research in our last common work, entitled: Nuclear Structure in Nonmesonic 
Weak Decay of ff „ W Q. has also been 8 n gg ested b y Dnbravfco Therefore it is not 
difficult to figure out who was the bird and who was the fish in our teamwork. 




FIG. 1: Escher's engraving where the entanglement between birds and fishes pictures my joint work with 
Dubravko. 
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Among several topics on weak interactions that I have tackled with Dubravko, I will limit 
the present discussion to a few features of the /3/3-decay, from which we can learn about the 
neutrino physics, provided we know how to deal with the nuclear structure. This is a second- 
order weak process whose electromagnetic analogies are the atomic Raman scattering and 
the nuclear 77-decay 

In nature there are about 50 nuclear systems where the single /3-decay is energetically 
forbidden, and therefore the /3/?-decay turns out to be the only possible mode of disintegra- 
tion. It is the nuclear pairing force which causes such an " anomaly" , by making the mass of 
the odd-odd isobar, (JV — 1, Z + 1), to be greater than the masses of its even-even neighbors, 
(N,Z) and (JV, —2, Z + 2). 

The modes by which the /3/3-decay can take place are connected with the neutrino (y)- 
antineutrino (v) distinction. If v and v are defined by the transitions: 

n —* p + e~ + v (1.1) 
v + n — > p + e~, 

the two-neutrino double beta (2v(3(3) decay (JV, Z) (N — 2, Z + 2) can occur by two 
successive (3 decays: 

(N,Z) ^ (N-l,Z+l) + e- + P 

^ (JV - 2, Z + 2) + 2e~ + 2z>, (1.2) 

passing through the intermediate virtual states of the (JV — 1, Z + 1) nucleus. 

However, neutrino is the only fermion that lacks a conserved additive quantum number 
to differentiate between v and v. Thus, it is possible for the neutrino to be a Majorana 
particle [v = v), i.e. equal to its own antiparticle a la 7r° . In this case the neutrinoless (3(3 
(OvfiP) decay is also allowed: 

(N,Z) ^ (N-l,Z + l) + e- + i> = {N-l,Z + i) + e- + v, 



0- 



(JV -2,Z + 2) + 2e~. (1.3) 



In absence of helicity suppression (as would be natural before the parity violation has been 
observed) this Ov mode is favoured over the 2v mode by a phase-space factor of 10 7 — 10 9 : 
T 0u rsj (10 13 - 10 15 ) yr, while T 2u ~ (10 20 - 10 24 ) yr. 



1 A Dirac particle can be viewed as a combination of two Majorana particles with equal masses and opposite 
CP properties, in which case their contributions to the §v(3fi decay cancel. 
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By the early 1950's several searches for the (3(3 decay were performed, inferring that 
T2v+ov ~ 10 17 yr. This pointed towards v ^ v, and prompted the introduction of the lepton 
number L to distinguish v from v: L = +1 was attributed to e~ and v and L = — 1 to e + 
and v. The assumption on conserving the additive lepton number allows the 2v(3(3 decay 
but forbids the 0^(3(3 one, for which AL = 2. 

But in 1957, with the discovery of parity non-conservation in weak interactions, the 
question of the Majorana/Dirac character for the neutrino is raised again. That is, the 
neutrino was found to be left-handed (LH) and the antineutrino right-handed (RH), and in 
place of (jl.2j) one now has: 

n -> p + e" + vrh, (1.4) 
i^lh + n — > p + e~. 

Consequently the second process in (jl.3|) is forbidden because the right-handed neutrino, 
emitted in the first step, has the wrong helicity to be reabsorbed in the second step. For 
massless neutrinos, as was believed they were, there is no mixture of vlh in vrh-, and the 
0i//?/9-decay cannot go through, regardless of Dirac or Majorana nature of the neutrino. This 
event discouraged experimental work for a long time. 

However, with the development of modern gauge theories the state of affairs began to 
change, and over the past decade the interest in the /3/3-decay sprang up again. There 
are many reasons for that. The most important one is the fact that the 0^/9/9-decay plays 
a decisive role in shaping the ultimate theory in any new physics beyond the standard 
SU(2)l x U(l) gauge model of electroweak interactions. Moreover, no solid theoretical 
principle prevents neutrinos from having mass, while the most attractive extensions of the 
standard model require neutrinos to be massive. Nor does theory predict the scale of neutrino 
masses any better than it can fix the masses of quarks and charged leptons. 

Yet, once the neutrino becomes massive, the helicity is not a good quantum number any 
more. Then, if the neutrino is in addition a Majorana particle with an effective mass (m^), 
the mixture of ulh in vrh is proportional to (m u ) / 'E u , and Ou (3 /3-decay is allowed 2 . This 
fact inspired experimental searches in many nuclei, not only for the Ou /3/3-decay but also 
for the 2z/-decay, since these two modes of disintegration are related through the the nuclear 

2 For simplicity we assume that weak interactions with right-handed currents do not play an essential role 
in the neutrinoless decay. 
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structure effects. In fact, their half-lives cast in the form: 



T,- 1 = G 2 uMl, T -J = g uMl(m u ) 2 , (1.5) 

where Q's are geometrical phase space factors, and A4's are nuclear matrix elements 
(NME's). M.iv and M.q v present many similar features to the extent that it can be stated 
that we shall not understand the Oz//3/3-decay unless we understand the 2i//3/3-decay. 

I will limit the discussion here to the 2v mode, but the whole presentation that follows can 
be straightforwardly applied to the Oz/ mode. With Dubravko we developed a full formulation 
of the O^/3/3-decay, based on the Fourier-Bessel expansion of the weak Hamiltonian, expressly 
adapted for nuclear structure calculations M, lal. We have also worked together on the 

fin 

"charged majoron models" |5|,|6|], so called because the majoron carries the unbroken £7(1) 
charge of lepton number. These models are probably the only ones that have a chance of 
producing the neutrinoless /3/3-decay that includes the emission of a massless majoron at a 
rate which could be observed in the present generation of experiments. 

At present we have at our disposal beautiful data on several 2i//3/3-decays, namely in the 
following nuclei: 48 Ca, 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te, 136 Xe, 150 Nd, 238 U and 
244 Pu. The 2z//5/5-decay turned out to be one of the slowest processes observed so far in 
nature and offers a unique opportunity for testing the nuclear physics techniques for half- 
lives ^ 10 20 yr. Disappointingly, as yet, after many years of heroic efforts of many physicists, 
no evidence of the nonstandard /9/?-decay has appeared. A survey of experimental results is 



given in the review article by Ejiri 



All we know is that the massive neutrinos can seesaw. In fact, the major advances in 
neutrino massiveness have been based so far on the compelling evidence that favors the 
neutrino oscillations. They first came from atmospheric neutrino flux measurements at 
SuperKamiokande (SK]_ lj| , and the solar neutrino shortages at SAGE, Gallex, GNO , 
Kamiokande and SK |l8|. However, only recently the solar SNO experiment [19(, jointly 
<!amLAND (KL) experiment 



with the reactor 
K2K experiment 



2Cj, and the first long baseline accelerator 



21], firmly fixed the neutrino oscillations. 
That is, we now know the neutrino mix, and we have initial values for their mixing matrix 
elements. We know the number of light active neutrino species and the differences between 
the squares of their masses. But we still don't know two crucial features of the neutrino 
physics, which are: a) the absolute mass scale, and b) whether the neutrino is a Majorana 
or a Dirac particle. Only the 0^/3/?-decay can provide this information, and this fact has 
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motivated the undertaking of very attractive next generation experiments for many different 
isotopes, including 48 Ca, 76 Ge, 100 Mo, 116 Cd, 130 Te, 136 Xe, 150 Nd, and 160 Gd 

The interest in neutrinos goes beyond the study of their intrinsic properties, and extends 
to a variety of topics in astro-nuclear physics, such as the understanding of the energy 
production in our sun, the synthesis of heavy elements during the r-process, the influence 
of neutrinos on the dynamics of a core-collapse supernova explosion and the cooling of a 
proto-neutronstar. The neutrino physics even appears in cosmological questions such as the 
role of neutrinos in the matter-antimatter asymmetry in the universe. 

The outline of this paper is as follows: In Section II we list a few general features of 
the NME's, which are necessary for understanding what follows. We discuss in Section III 
the general behavior of the /3/3-decay amplitude within the charge-exchange QRPA, and we 
record the matching refinements as well. We develop in Section IV a simple nuclear model 
for the /3/3-decay, based on the well-known Quasiparticle Tamm-Dancoff Approximation 
(QTDA). In Section V we do the analysis of the 48 Ca^-> 48 Ti 2z/-decay in the particle- hole 
limit of this model. A few final comments and remarks are pointed out in Section VI. 



II. 2upp MATRIX ELEMENT 



Independently of the nuclear model used, and when only allowed transitions are consid- 
ered, the 2v(3(3 matrix element for the |0t) final state reads 

'<o;i|Ori|A+>(A+||Of||0+>' 



M 2v (f) = EH A E 

A=0,1 a 



V 



MUf)+M%(f) (2.1) 



where the summation goes over all intermediate virtual states |A+), 



Of = (2A + l)- 1 / 2 E(p||OA||n) (4cn) with 

pn 



=1 for F 

01 = a for GT 



(2.2) 



are the Fermi (F) and Gamow- Teller (GT) operators for /3~-decay, and (c) are the particle 
creation (annihilation) operators. The corresponding /5 + -decay operators are = (O^ ) , 
and 



*D\+ f = E\+ — 



E + E + E + - E 

J TP TP J 



(2.3) 



2 2 

is the energy denominator. E and E + are, respectively, the energy of the initial state |0 + ) 
and of the final states |0j). 

6 



Contributions from the first-forbidden operators, which appear in the multipole expansion 
of the weak Hamiltonian, as well as those from the weak-magnetism term and other second 
order corrections on the allowed 2z//3/3-decay are not relevant for the present work and will 
not been tackled here. In recent years we have examined all of them rather thoroughly 

In nuclear physics the isospin symmetry is conserved to a great extent, while the Wigner 
SU(4) symmetry is not. Because of this, the amplitude M.^v * s often neglected in the 
literature. Nevertheless, one should keep in mind that, while the mean field strongly breaks 
the isospin symmetry, the residual force restores it almost fully. Therefore, although in many 
cases the final value of M.^ v is small, it is recommendable to keep track of this NME during 
calculation so as to test the consistence of the nuclear model, as well as to fix its coupling 
constants. 

The energy distributions of the transition strengths |(AJ||0f ||0 + )| 2 link the sin gle (3 ± - 
decays to the single charge-exchange reactions, such as (p,n), (n,p), etc jlfl I25I I26II27I l^sj 
Q. The total /3 ± strengths 

sf = (2A + l)- 1 £i<Aj||<ijO+}| 2 , (2.4) 

a 

can be expressed in the form 

sf = (0 + |Of -Of |0+> = (-) A (2A + l)- 1 (0 + |[Of of ] |0+>, (2.5) 

when |A+) is a complete set of excited states that can be reached by acting with Of on the 
initial state |0 + ). It follows at once that 

Si = Sf - Sf = R A (2A + lrWtOf , Of ]o|0 + > =N-Z, (2.6) 

which is the well-known single-charge-exchange sum rule, also called Ikeda sum rule (ISR) 
(25!, for both the F and the GT transitions. 

Similarly, the /3/3-decays are closely related to the double-charge-exchange reactions, and 
to the spectral distribution of their strengths, 

Sf^f) = (2A+l)- 1 |E(0;i|Of ± ||A+)(A+||of |0+>| 2 , (2.7) 

a 

over the final states |0t). The total strengths are defined as: 

Sf" = E^(/) = (2A+l)- 1 EK0?|Of 'Of l° + )| 2 (2-8) 



and can be rewritten in the form 

Sf ± = (2A + I)" 1 (0+| Of • Of of • Of |0 + ), (2.9) 
The double-charge-exchange sum rules (DSR) are: 

gfifi = gPfT _ gPf* = ^ + l)-l <0 +| [O f . f , f . Of ]|0 + >, (2.10) 

which when evaluated give jjic], | : 

5^ = = 2{N - Z){N -Z-l), (2.11) 

and 

S^ T = Sf* = 2(N-Z)(N-Z-1 + 2Sf) - ~C, (2.12) 



where C is a relatively small quantity and is given by |3_1, (5)]. 

III. CHARGE-EXCHANGE QUASIPARTICLE RANDOM PHASE APPROXI- 
MATION AND BEYOND 

The (5(5 decays occur in medium-mass nuclei that are often far from closed shells, and, as 
a consequence, most of the recent attempts to evaluate M.2u and M.q v rely on the neutron- 
proton QRPA, because this model is much simpler computationally than the shell model 
(SM). Note that the kind of correlations that these two methods include are not the same. 
The QRPA deals with a large fraction of nucleons in a large single-particle space, but within 
a modest configuration space. The shell model, by contrast, deals with a small fraction of 
nucleons in a limited single-particle space, but allows them to correlate in arbitrary ways 
within a large configuration space. 

The charge-exchange QRPA has been first formulated, and applied to the allowed (5- 
decay and to the collective GT resonance, by Halbleib and Sorensen in 1967 However, 
intensive implementations of QRPA to /3/3-decay began only about 20 years later when Vogel 



and Zirnbauer 



discovered that the ground state correlations (GSC) play an essential role 



in suppressing the 2v(5(5 rates. Soon afterwards, Civitarese, Faessler and Tomoda 



to the same conclusion. Almost simultaneously, Tomoda and Faessler |35[ , and Engel, Vogel 



34 1 arrived 



and Zirnbauer |36j revealed a similar though smaller effect on the 0v(5(5 decay. 

When applied to the /3/3-decay the following two steps are performed within the standard 
QRPA: 
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1. Two charge-exchange QRPA equations are solved for the intermediate 1 + states; one 
for the initial nucleus (N, Z) and one for the final nucleus (N — 2, Z + 2). The first 



one, 



A B 
B A 

is evaluated in the BCS vacuum 



X 
Y 



UJr, 



X 

-Y 



(3.1) 



|0 + ) = II( M P + v P C t C p) IlK + v n4i c V)\), 



(3.2) 



where |) stands for the particle vacuum. The Eq. (|3.1|) describes simultaneously four 
nuclei: (N-1,Z-1), (N + 1,Z-1), (N-1,Z+1) and (JV+ 1, Z + 1). The matrix 
elements of the operators Of are: 



pn 

(UII<l|o + ) = E[ l0 -(^;i)^n ; i+ + i+(^;i)^ ;1 H 



(3.3) 



were 



A° + (pn;X) = ttpVft(p||0 A ||n), 
A^(pn;X) = v p u n (p\\Ox\\n), 

are the unperturbed strengths. The ISR reads 

^T = ^E[l(Ull^nio + )l 2 -KUII<llo- 



(3.4) 



N-Z. 



(3.5) 



In the same way the second QRPA does not deal with the (N — 1,Z+1) nucleus only, 
but entangles as well the isotopes (N — 1, Z + 3), (N — 3, Z + 1) and (N — 3, Z + 3). 
The corresponding ISR is 

%T = \i: [I CC 1 1 Of" 1 10 + ) | 2 - | (i: \\Ot l|0 + )| 2 ] =N-Z = N-Z-A 1 (3.6) 

where the barred kets (|0 + ), |1„)) indicate that the quasiparticles are defined with 
respect to the final nucleus. 

2. The equation (j2.1|) is substituted by one of the following two ansatz: 
• Method Ml (proposed by Vogel and Zirnbauer 0,Q]): 



9 



= -E 



(i+ii<iio + )(uii^ iiq + ) , (Tiii^f 1 ir)(Tiiigrno") 



^ + iin+\/T+ll^i/3-|ln^ 



Method M2 (introduced by Civitarese, Faessler and Tomoda 

(Tiii<no + )(i:ii c t)(uiioriio + ) 



, (3.7) 



M 



2p 



2E 



U-.+ + uj T + 



(3.8) 



where the overlap in ()3.4|) is evaluated as: 



pn 



(3.9) 



Note that the last two equations for M.2vi ()3.7|) and ()3.8|) . cannot be derived mathe- 
matically, and that they are just recipes which make the applications of the QRPA to the 
/3/3-decay possible. Moreover, in many applications of the method Ml, the energy denomina- 
tor (co>j+ +ujj+ f )/2 in ()3.4|) is simply taken to be equal to u) x +. The GSC in (l+||Of ||0 + ) and 
(1^/| \Oi + 1 10 + ) are also different, matching, respectively, transitions (N, Z) (N+l, Z — 1) 
and (JV-2,Z + 2) ^ (N-3,Z + 3). 

When compared to the SM, the QRPA presents the following drawbacks: 

I) There is ambiguity in treating the intermediate states, and further developments must 
be made to match the excited states of the odd-odd nuclei based on different ground 
states of the initial and final even-even nuclei, as in (J3.7)) and (J3.8I) . 

II) The QRPA /3/3-decay amplitudes are extremely sensitive to the PN interaction in the 



particle-partic 
S = forces 



yp) channel, or more precisely to the ratio between the S = 1 and 



33 



34. 



35 



37 



38, 



39 



40 



41 



42 



43|. 



Even worse, the model 



collapses as a whole in the physical region of t or g 



pp 



III) The QRPA is only capable to account for the /3/3-decay into the ground state |0^) 
of the final nucleus, while to deal with final excited states it is necessary to recur to 
another nuclear model. Usually, one solves an extra charge-conserving QRPA equation 
of motion for the 2 + -excitations on the final BCS vacuum, and one assumes the final 
states |2+) and |0^~) to be the one-phonon quadrupole vibration, and a member of the 
two-phonon quadrupole vibrational triplet, respectively 



45 



46 



47 



48, 



49]. 



For a 5 force this ratio is t = v\ v jv 1 ^ . Usually is introduced in an ad hoc parameter, denoted by g : 



pp 
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IV) In the QRPA we cannot evaluate the energy distributions of the double charge- 
exchange transition strengths, given by ([2.7)1 . 

There have been many tries to avoid the collapse in the QRPA, so its sensitivity to the 
phenomenological par ameters would become more realistic. A first step in this direction has 
been done in Ref. 37 1, shortly after the finding of Vogel and Zirnbauer j^, by working out 
a two-vacua QRPA (TVQRPA) specially tailored for the /9/9-decay. We called it so because 
one solves the QRPA equation of motion on the quasiparticle vacuum 

|o + ) = IIK + vM) Il(«n + ^44)1), (3.io) 

p n 

which involves the vacua of both the initial (it, v) and final (u, v) nuclei. 

The forward (A) and the backward (B) going matrix elements are given by 

A(pn,p'n';J) = (e~ +e^)5 pn>p , n , 

+ s/PpPnPp'Pn' [(u p v n Up>v n > + v p u n Vp>u n >)F(pn, p'n'; J) 

+ (u p u n u p 'U n i + v p v n v n >u n r)G(pn,p'n'; J)], (3.11) 

and 

B(jm,p'n'; J) = y/p p p n p P >Pn> [(v p u n u p >v n > + u p v n v p >u n >)F(pn,p'n'; J) 

- (u p u n v p >v n > + v p v n u p >v n >)G(pn,p'n'; J)} . (3-12) 



The quasiparticle energies are: 



2u p v p pp 2u n v n p n 

with the pairing gaps given by 



(3.13) 



A ™ = ~\ \ 2 JW/ + l u n>Vn>Pn>G{nn, n'ri; 0). (3.14) 

The factors p p and p n are defined as p~^_ = u p + v p , p" 1 = v? n + u^, while the remaining 
notations have standard meanings j^, 13 ■ 
The M.iv amplitude reads 

M ^ = E MW) (315) 
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FIG. 2: Graphical representation of the numerator in (|2.1() within the QRPA. The first and second terms 
match, respectively, to (0 + and (A + ||0^ ||0 + ). The residual interaction brings about the ground- 

state correlations in the vertex V. 

where 

<^H<?niO + > =E[At(wl)X pn;li +A _(wl)W , 

pn 

<^ll<l|0+> = £ [^(^1)^+^(^1)^-1+] > (3-16) 

pn 

are, respectively, the perturbed (3~ and /3 + matrix elements, and 

A° + (pn;X) = y/PpPnU p v n {p\\Ox\\n), 

A°_(pn;\) = y/p p PnVpU n (p\\Ox\\n), (3.17) 
are the unperturbed ones. The ISR is now 

%r = \ l2[\( l t\\Or\\0 + )\ 2 - |<1+||<||0 + >| 2 ] =N-Z-2. (3.18) 



Now and then the gap equations are solved for the intermediate nucleus in which case the 
ISR gives exactly - Z - 2 flfl0. 

In Figure 2 we exhibit the graphical representation of different terms within the numer- 
ator in (|3.13j) . The TVQRPA involves only the virtual (N — 1,Z + l)-states, being the 
backward-going contributions for the j3~ transitions the forward-going contributions for the 
j3 + transitions, and vice versa. 
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Based on a numerical comparison, we have found that the behavior of 2u transition 
amplitude in the TVQRPA, with regard to the pp coupling, is quite similar to that previously 
noticed in the QRPA j^. Since no uncertainties are present when intermediate states are 
treated in the TVQRPA, this result can, in some ways, be interpreted as a justification for 
the averaging procedure performed in Eqs. ()3.7|) and ([3.8)1 . 

It could also be worth noting that in the TVQRPA we can express the amplitude M.% v 
in the form: 



M 



2u 




AO 
-A 



(3.19) 



and can therefore calculate the transition probability without first solving the QRPA equa- 
tion. 

We found the last equation to be especially useful for discussing the /3/9-decay rates within 
the single-mode- model (SMM) Q,HQ, 

which deals with only one intermediate state for 
each J n and is the simplest version of the QRPA. Within the SMM one can express the 
moment A4 2 u a la Alaga 



M 2v = g e 2 f jMt, 



(3.20) 



i.e. as the unperturbed BCS matrix element 



(p\\Oi\\ny 



M\ v = u p v n u n vpp p p n - n 



(3.21) 



multiplied by an effective charge: 



9-iv 



Co>iH 



1 + 



G(pn,pn; V 



(3.22) 



where 



u = — - 



4 



G(pp,pp; + ) + G(nn, nn; 0" 



(3.23) 



is the unperturbed pairing energy between protons and neutrons. g^J comes from the QRPA 
correlations, or more precisely from the interference between the forward and backward going 
contributions, which add coherently in the pp channel and totally out of phase in the ph 
channel. 

Still more, as the QRPA energy ui+ is in essence a linear function of G(1 + )/uj°, we can 
state that, because of the GSC, the effective QRPA charge g 2 J is mainly a bilinear function 
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of G(l + )/u°, or equivalently of t |40j,|4 



i.e. 



&' ~ (3-24) 

That is, A^2f passes through zero at t = to where G(l + ) = —uj°, and has a pole at t = to 
where uji+ = 0. 

The fact that in many situations the SMM reproduces quite well the complete QRPA 
calculations made us suspect that the behavior of M. 2v with regard to t could always follow a 
(1, 1)-Pade approximant of the form ()3.21j) . Thus, we have suggested that, independently of 
the nuclear Hamiltonian and/or the configuratio n sp ace employed in the QRPA calculation, 



the 2i> amplitude should unavoidably behave as 



4( 



M 2u = M 2u (t = 0)^—^. (3.25) 

l — t/tx 

At variance with the bare BCS value M 2u , given by (|3.20|) . the matrix element M. 2v it = 0) 
also contains the ph-like correlations. We have tested the relation ()3.24|) only for a zero- 
range force 40], but as far as I know there is no QRPA calculation in the literature that 
could be in conflict with this result . 

One can also guess that the breakdown in the QRPA comes from the violation of the 
particle-number symmetry, caused by the BCS approximation. Therefore with the hope of 
getting out of this inconvenience we have worked out a particle number projected QRPA 
(PQRPA) for the charge-exchange excitations starting from the time-dependent variational 
principle, j^. However, after performing numerical calculations for the 2z//3/5-decay in 76 Ge, 
we found that in this model the Ai 2u amplitude continues to behave roughly as in the plain 
QRPA. Said in another way, the number projection procedure is unable to avoid the collapse. 



The variation of the QRPA that 



renormalized QRPA (RQRPA) 



50 



ras received major attention lately is the so-called 



51 



52 



53, 



54 



55 



56 



571 ] . The new ingredient brought 



up by this model is the effect of the GSC in the QRPA equation of motion itself. The 
important outcome of this is that the QRPA collapse does not develop anymore in the 
physical region of the j>j>-strength parameter. Yet, this new procedure to incorporate the 
GSC tones down only slightly the strong dependence of the 2v/3/3 transition amplitude on 
this parameter. 



1 When the renormalization coupling constant g pp is used j34j a similar expression to the Eq. I|3.22|) is valid 
(with gPP's for i's). 
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On the other hand we soon found |51| that the price one has to pay in the RQRPA 
to avoid the collapse was the non-conservation of the ISR, given by the Eq. ()2.6j) . This 
violation is about 20 — 30% and we cannot get away from it. It comes from the fact that the 
scattering part of the GT operator, when acting on the RQRPA ground state, creates states 
that are not contained in the model space. These terms have recently been considered in 
the framework of the "Fully Renormalized QRPA" (FR-QRPA) 58], whereby the ISR was 
successfully restored. Yet, within the FR-QRPA the 2uf5f5 amplitude behaves similarly as in 
the ordinary QRPA. Namely, in this model M.2v passes through zero and develops a pole for 
values of the pp-strength parameter which are only slightly higher than those in the QRPA 
model. 

Let us also remember that only the self-consistent QRPA (SCQRPA) theory 



59 



60 



61 



62( incorporates fully the GSC, leading simultaneously to a coupling of the single-particle 



field to the QRPA excitations. We performed [61[ a detailed comparison of the properties 
of the QRPA, RQRPA and SCQRPA equations in the 0(5) model for the F excitations, 
inferring that: i) before the QRPA collapses all three approaches reproduce well correct 
results, ii) near the transition point only the SCQRPA values are close to the exact ones, 
and iii) beyond that point both the SCQRPA and RQRPA yield values different from the 
exact ones, but the former are somewhat better. One can suspect that in realistic cases this 
condition prevails and, to some extent, even for the GT transitions. Such a possibility has 
not been explored so far. It should also be mentioned that what some authors j6jj] refer to 
as self-consistent QRPA is just the RQRPA with the introduction of some minor changes, 
given by Eqs. (34), (40) and (45) in Ref. 4 

We arrive therefore to the conclusion that not one of the amendments of the QRPA, 
proposed so far to rescue this nuclear model, was able to change qualitatively the behavior 
of the amplitude M.2v, given by ()3.22|) . unless we agree to tolerate the violation of the ISR, 
which could be extremely dangerous as we have no control on how it affects the definite 
value of Ai2w More, neither the RQRPA nor the SCQRPA is able to evade the other three 
unfavorable QRPA outcomes that we have pointed out at the beginning of this section. Thus, 
within this QRPA scenario, instead of introducing further improvements and variations into 
the QRPA equation of motion, it would perhaps be a good idea "to shuffle the cards and 
deal again" , i. e. to try to work out a different quenching mechanism for the NME's. In the 



See also comments with regard to this in Ref. |64j . 
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next two sections we explore that possibility. 



IV. QUASIPARTICLE TAMM-DANCOFF APPROXIMATION (QTDA) 



Here we sketch a simple nuclear model for evaluating the (3(3 decay rates, based on the 
well-known QTDA 65J, |66j. It will become clear immediately that the main difference in 
comparison to the QRPA comes from how one describes the final (N — 2, Z + 2) nucleus. 

Same as in the QRPA, we conveniently express the total Hamiltonian as 



H = H v + H n + H pn + H pp + H nn = H + H,, 



(4.1) 



where H p and H n are the effective proton and neutron single-quasiparticle Hamiltonians 
(with eigenvalues e p and e n ), while H pn , H pp , and H nn are the matching effective two- 
quasiparticle interaction Hamiltonians among the valence quasiparticles. 

We assume both i) that the initial state is the BCS vacuum in the (N, Z) nucleus, and 
ii) that the intermediate and final nuclear states involved in the /?/3-decay are, respectively, 
two and four quasiparticle excitations on this vacuum. That is: 

• initial state: |0 + ) = \BCS), 



intermediate states: 



with 



\K) = Y, X pn;\+\P n -> x+ )i 



pn 



final states: 



with 



\pn;X + ) = [alaU x+ \BCS), 



l°/>= E ^ip 2 nm 2 J;0+ \PiP2, n\Tlz\ J) , 
PiP2n\n 2 J 



|piP2,nin 2 ; J) = A^(p 1 p 2 )iV(nin2){[4 1 aJ 2 ]j[4 1 aJ l2 ]j}o|-BC i S) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



and 



N(ab) = (l + 5 ab y 1 / 2 . (4.6) 
Here (a) is the quasiparticle creation (annihilation) operator relative to the BCS 



vacuo. 



16 



FIG. 3: Graphical representation of M.2v in the QTDA. The first and the second vertex match with matrix 
elements l|4.7|l and Q4.8p . respectively, while the third vertex represents the residual interaction in the final 
state. 

One can read the matrix elements of H pn between the intermediate states \pn; A + ) and the 
final states \p1P2, nin 2 ; J), respectively from |66|, (4.9) and (4.5)], and those of H pp and 
H nn between the same final states from jfisl (2.11)]. We show explicit results only for the 
one-body matrix elements that appear in (|2.1|) . They are: 

<A+||Of l|0 + > = £A^;A)X pn;A +, (4.7) 



and 



(0;i|OTl|A+> = -E^ P n;A+ E ^ piP2niI12 j;o+^(^2)iV(n 1 n 2 )P(p 1 p 2 J)P(n 1 n 2 J) 

pn PiP2niri2J 

x v ^7TT(-r +P2+J+A I Pl Hl X \A° + (p 1 n 1 ;\)5 P2P 6 n2n , (4.8) 



p n J 



where 



p(p lP2 j) = 1 - (-r +p " +j p( Pl <- P2 ), (4.9) 

is the well known permutation operator. The energies in the denominator (|2.1|) are 

E \+ = E o + + \> ~ K> 

E Q + = E + lu + + 2\ p - 2A n , (4.10) 
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where uj x + and u + are the eigenvalues of the Hamiltonian (j4.1|) for intermediate states |A+) 
and final states |0t), respectively, and A p and A n are the chemical potentials. Therefore 

The QTDA has the correct particle- hole (shell model) limit: v p — > 0, v n — > 1, and therefore 
one can straightforwardly apply this model to single- and double-closed shell nuclei. In the 
next section we discuss an example. 



V. 2^/3/3-DECAY 4S Ca 



is 



Ti 



There are two recent experimental results for the 2z//3/3-decay half-life to the Of state 
that nicely agree with each other. They are: 



Ref. 
Ref. 



68j : 



T 2v = (4.3+^[stat] ±1.4[syst]) x 10 19 yr, 



T, 



2v 



A 0+3.3 



x 10 19 yr 



which, from (jl.5|) and the kinematical factor [4 

Q 2v = 42.3 x KT 19 

yield: 

\M 2 „\ = 



yr(MeV)' 



(5.1) 



(5.2) 



Ref. 
Ref. 



0.07418;^) [MeV]- 1 



\M 2u \ 



0.075±8:8ii) [MeV]- 1 . 



(5.3) 



One should keep in mind that the experimental values for |A^2v| i n (|5.3p depend, through 
the value of Q 2u , on the value used for the effective axial- vector coupling constant qa- In the 
present work we use gA = 1- Clearly for the bare value, gA = 1.26, the phenomenological 
NME's decrease by factor (1.26) 2 . 

n 

There is also a very interesting high- resolution charge-exchange reaction experiment [69J , 
where M. 2v for the ground state in 48 Ti, was built from energy spectra of the 48 Ca (p, 
n) 48 Sc and 48 Ti(d, 2 He) 48 Sc reactions, by converting the (p,n) and (d, 2 He) cross sections 
into moments (0i"||Of an d (lJ||Oi ||0 + ) which contribute in (|2.1|) . In performing the 

summation over a five experimentally observed states below 5 MeV have been considered, 
under the assumption that all matrix elements are positive. In this way Rakers et al. 
get: 



Ref. 



\M 



2v\E<bMeV 



(0.0740 ±0.0150) [MeV]" 1 . 



(5.4) 



18 



One should 



and that 



although (jO □ 



50 




E (MeV) 



FIG. 4: Graphical representation of energy dependence of the double GT strengths Sqj, , measured from 
the 48 Ca ground state. The coupling constant values in the perturbed calculations were v s = 35, Vt — 65 
for case (a) , and v s = 40, vt — 60 for case (b). 



Previously, jointly with Dubravko, we had already studied the (3(3 decay in 48 Ca, but only 
for the single particle state I/7/2 Here we consider the complete pf particle-hole space. 
Of course, this configuration space is still stron gly limited when compared to those of the 



SM within the pf single-particle space iTOt l7ll . l72j ]. Thus, obviously, a SM yields a more 



realistic description of the /3/3-decay [6JJ,|68[ and of the related charge-exchange reactions 
than the present model. Our only purpose here is to explain the way the QTDA quenching- 
mechanism works. Choosing suitable effective single-particle energies (SPE) is, as always, a 
delicate issue in nuclear structure calculations. We will take them to be: 



€ jn e j n + ^jr, 



where e- n are experimental SPE for Ca, extracted from Fig. 2 in Ref. 



(5.5) 



namely: 



z h/2 



' J 5/2 



6.5, el 



2-1, and e° i/2 



4.1, in units of MeV. The self-energies /jLj n and 
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. 3- 



Ac is the Coulomb 



Hj p account for 8 neutrons in the f 7 / 2 shell, and are defined in Ref. 
displacement energy in A8 Ca. With SPE chosen in this way the energy of the isobaric analog 
state (IAS) is always equal to Ac, as it should be. 
We will use the 5-force 

V = -4n{v s P s + v t P t )5{r) (5.6) 

both for evaluating the self-energies in ()5.5|) and for calculating the residual interaction. 
Besides, for the purpose of simplifying the discussion, we will assume the coupling strengths 
v s and v t to be equal for identical and for different particles. 

Let us first comment the unperturbed results, i.e. when v s — Vt — 0. In this case 
all double F /3/3-strength, Sp = Sp 13 = 124, concentrates in the lowest lying degenerate 
states \f 7/2 , f 7/2 ] J = 0,2,4,6), in parts of S p /~ {J) = 4(2 J + 1). These are the double 
IAS's (DIAS's) and are at energy 2Ac = 12.78 MeV, measured from the 48 Ca ground-state, 
which was taken to be E = 0. Contrarily, we found in these states only 12% of the total 
double GT intensity, S^f T = SgfT = 125.71, being equal to: 2.20, 7.22, 2.64 and 3.18 for 
J = 0,2,4 and 6, respectively. The major part of the double GT strength concentrates at 
levels |/7 /2 , f 7 / 2 f 5/2; J = 2, 4, 6) (55%) and \ ff /2 , /| /2 ; J = 0, 2, 4) (33%), which, as seen from 
Figure 4 lie, respectively, at energies 2 A c + e/ 5 2 , and 2A C + 2e} 5 2 - For the lowest final 0+ 
state the unperturbed denominators ()4.10|) are all null, which makes both NME's, M. 2v and 
M-2v 1 t° become 00. The scene changes radically when the residual interaction is switched 
on. First, as we show in Figure 5, there is a great variety of physically sound values for v s 
and Vt that allow the model to account for the phenomenological NME's f!5.3|) and (|5.4j) . 
In another words, same as the QRPA, the QTDA is capable of restoring the Wigner SU(4) 
symmetry, quenching in this way the NME's. This is precisely the mechanism we have been 
searching for. 

Just for the sake of illustration we use here two sets of coupling constants: (a) v s = 35, 
v t = 65, and (b) v s = 40, v t = 60, which reproduce reasonably well t 
shown in Figure 6, and are close to values used in our previous works 
PN ph channel (vf = 27 and vf = 64). Note first that, as (f 7/2 , / 5/2 ; l + \H pn \f 7/2 , f 7/2 ; 1+) 
(v s + Vt), the wave function for the GTR is the same in both cases: 



re energy levels of 48 Sc, 
Q,Q,0 within the 



I GTR) = 0.972|/ 7/2 , / 5/2 ; 1 + ) - 0.237|/ 7/2 , f 7/2 ; 1 + ). (5.7) 

The residual interaction acts differently on double F and on GT transition strengths. In the 
first case, all intensity remains concentrated at the energy 2A C , independently of the values 
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FIG. 5: Graphical representation of the (3(3 amplitudes M§J in the QTDA as a function of coupling 
constants v s and v t . 

of v s and v t , but now there is only one DIAS. Its wave function is a coherent superposition 
of the states \fy 2 , fy 2 ; J = 0,2,4,6), i.e. 



|dias) ^VeOT 



1/7/2) /7/2? J)- 



(5.8) 



In the second case the double GT resonance (DGTR), placed at an energy that is only slightly 
lower than 2Ac + e} , carries a large part of the strength. That is: S^ T (DGTR)= 72.87 
in case (a) and = 72.61 in case (b). The DGTR wave function is: 



|DGTR) 



0.405 
0.371 



}\fr/2,f^0) + 



0.683 
0.683 



* 1/7/2) f5/2i 2) + 



0.505 
0.533 



I/7/2, / 5 2 /2 ; 4) + 



(5.9) 



The remaining GR strength mainly comes from the states \ fy 2 , ^2/5/2; J — 2, 4, 6), and in 
both cases concentrates almost fully at the energy of ~ 10 — 11 MeV. The unperturbed and 
perturbed double GT strengths are confronted in Figure 4. The calculated wave function 
for the ground state in 48 Ti is 



|0f 



0.922 
0.922 



M/7%/7 2 / 2 ;0> - 



0.250 
0.226 



I/7/2? fy2'i 2) - 



0.095 
0.093 



M/7/2,/7 2 / 2 ;4) 
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1.5 - 



> 
0) 



LU 

< 



0.5 




- 1.5 



0.5 



FIG. 6: Comparison between the experimental and theoretical energy levels of 48 Sc for the cases (a) v s = 35, 
v t = 65, and (b) v s = 40, v t — 60. We also display the shell-model results from Ref. |71|. 



0.189 _ _ 0.119 , 

+ l/7 2 /2 ; /5/2;0)+ i/7/2,P3 /2 ;o} + ---. (5.10) 

0.211 J 0.138 1 

It appears at the energy E Q + = —4.41 MeV in case (a) and at —3.95 MeV in case (b), which 



67 



68, 



compares favorably with the measured value Qpp- = 4.27 MeV 

The way in which Mfj(0f) are constructed from individual transitions to the interme- 
diate states |1+) is shown in Table I. We see that the GTR contributes either destructively 
- as in case (a), or its contribution is negligibly small - as in case (b). The first result is 
consistent with the SM calculation j^. The interference effect is still more pronounced on 
the GT strength going to the ground state in 48 Ti. In fact, for all practical purposes it turns 
out to be null in case (a): S^ T ~{0f) = (0.198 - 0.176) 2 /3 = 0.000, and extremely small 
in case (b): Sqj, (®f) = 0.036. Note, however, that the quenching mainly comes from out 
of phase contributions among the seniority-zero and seniority-four components in the wave 
function (|5.1()jl . which is responsible for relatively small values of the moments (0^ \ \Of \ |1^) 
shown in Table I. 

The calculated GT, F and total 2i>f3f3 matrix elements are shown in Table II. One sees 
that the F moments are quite significant and that they cannot be neglected as is usually 
done. In point of fact, in several calculations, where sizable contribution to the total M.$ y 
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TABLE I: Decomposition of the total GT matrix elements. The denominators T),+ , are given in 
units of MeV, and the moments M 2 J m units of (MeV) -1 . 



a 



<l+||0f~||O+> (0+||On|l+) (l+||0ri|O+)-(O+||Of"||l+) V lU M c 2 



iGT 



2.240 
-4.357 
total 



2.240 
-4.357 
total 



case (a): 
0.088 
0.041 



case (b): 
0.142 
0.002 



v s = 35, vt = 65 
0.198 
-0.176 



Vs = 40, v t = 60 
0.318 
0.010 



3.40 0.058 
11.32 -0.016 
0.042 



4.01 0.079 
11.91 0.001 
0.080 



TABLE II: Results for GT, F and total 2v(3(3 matrix elements in units of (MeV) -1 . 



case 


Mg T 


Ml 


M 2v 


(a) 


0.043 


0.008 


0.035 


(b) 


0.080 


0.011 


0.069 



moment was found to come from the virtual + states, the F contribution to the M. 2v 
moment had been omitted. This is manifestly inconsistent. 5 

Similarly, the model wave function for the excited state | O2") in 48 Ti is basically composed 
of seniority-four basis states, i.e. 

0.430 1 _ _ 0.781 1 _ _ 0.373 1 _ 

|0+> = |/ 7 2 /2 , / 7 2 /2 ; 2) - |/ 7 2 /2 , / 7 2 /2 ; 4) + |/ 2 /2 , / 2 /2 ; 6) + • • • . 

0.504 J 0.760 J 0.307 J 

(5.11) 

Unlike for the ground state, in case (a) appears an appreciable amount of the GT strength, 
Sqt (O2) = 0.170, in the state \0 2 )- Nevertheless, we will not discuss here the /?/9-decay rate 

5 A more careful study should go beyond the allowed approximation considered here and include the higher 
order contributions in the weak Hamiltonian as well 8, 241. 
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for the exited + state because the model does not reproduce satisfactorily the excitation 
energy; namely, we get 5.72 MeV, while the experimental value is 3.00 MeV. Finally, we 
note that in case (b) Sgf ((#) = 0.038. 



VI. CONCLUDING REMARKS 



The QTDA does not suffer from the inconveniences that have been listed in Section III 
in relation to the QRPA. More specifically, the similarities and the dissimilarities between 
the two models are: 

1. While the QTDA contains two ll f3~- like" vertices (see Fig. 3), in the QRPA we always 
approximate one of them by a like" vertex (see Fig. 2). This statement is valid 
for all variations of the standard QRPA, such as the TVQRPA, PQRTA, RQRPA and 
SCQRPA. 

2. The QTDA moments (A+||0f ||0 + ), given by l)4.7J) . produce the F and GT resonances 

in the same manner as in the QRPA. Furthermore, as the backward-going QRPA 
contributions have rather little impact on the tl /3~- like" transition strength, given by 
the first equation in ([3.5)1 . both models yield very similar results. 

3. Similarly, moments (0+1 |Of | |A+), given by ()4.8|) . are strongly reduced by the residual 
interaction, the QRPA moments (A+||Cf ||0+), given by the second equation 
in ()3.5j) . restoring in this way the isospin SU(2) and Wigner SU(4) symmetries, broken 
initially by the mean field. We know that in QRPA this symmetry-reestablishment 
takes place through the cancellation effect between the forward and the backward 



going contributions. In fact, in several works |8j,|29|, 3,|43j| we have used the property 
of maximal restoration of the SU(4) symmetry to fix the value of the pp parameter t. 
Instead, in the QTDA the quenching comes from out of phase contributions among 
seniority-zero and seniority-four configurations in the wave function of the final state 

|0f>- 

4. In QTDA, unlike in QRPA, the intermediate states |A+) have a unique meaning. 

5. The QTDA never collapses, and therefore the amplitude M.iv does not have a pole 
anymore, as happens in ()3.24|) . So, one can expect that, as far as model parameters 
are concerned M.2v will behave more moderately in the QTDA than in the QRPA. 
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6. In the QTDA one obtains the transition amplitudes M.2v for the excited Oj states 
without any additional effort. One can also easily evaluate the decays to the ex- 
cited final 2^ states; it is sufficient to diagonalize the Hamiltonian (|4.1j) in the 
\{PiP2)Jpi { n i n 2) Jn\ 2 + ) basis, and to calculate the matching transition operator (J4.8|) . 

7. In the QTDA the energy distributions of the double GT transition strengths ()2.7|) 
can be evaluated directly and, the corresponding sum rule, given by (j2.12j) . could be 
occasionally violated to some extent. Contrarily, the Ikeda sum rule, given by (J2.6|) . 
is always fully conserved in this model. 

The remarks quoted above suggest that perhaps the QTDA might be a more "natural" 
and a more suitable nuclear structure framework for describing the /3/3-decay than the QRPA. 
In fact, one should mention again that the QRPA was originally formulated for the single 
^-decays j^, and only later adapted for the /3/3-decays via the Ml and M2 Q Q, S 
361 ] ansatz. The TVQRPA, which was mathematically tailored, specifically for the j3(3- 
decay, is free of these averaging procedures and has the correct BCS limit given by (|3.21|) . 
Nevertheless the latest model has received rather little attention in the literature. What's 
more, quite recently it has been claimed that the ove rlap (KID I) gives rise to an additional 
suppression mechanism for the NME's M.2u and M.q v |74|, |75J]. We fully disagree with such 
a view. 

We do not suggest that the QRPA should be substituted by the QTDA. We merely state 
that the use of these two nuclear models in a joint manner should very likely reduce the 
uncertainties in the evaluation of the NME's, which is at present one of the principal worries 
in the nuclear physics community, and which has engendered a great deal of activity in 



recent years [76 



77 



78,[7^. One cannot but highlight the Suhonen's article [3] where he 
argues that within the QRPA it is not possible to account simultaneously, i.e. with the 
same set of model parameters, for the simple and double /3-decays. 

In summary, in this work we have demonstrated that the simple version of the QTDA 
proposed here is able to account for the suppression of M.2vi which was the major merit 
of the QRPA. Moreover, we feel that this model comprises all essential nuclear structure 
ingredients that are needed for describing the /3/3-decay processes. Of course, whether this 
is totally or only partly true has still to be tested, and it might be convenient to consider 
some additional refinements in the future. Their incorporation, however, should not, in 
principle, involve serious difficulties. For instance, the number projection procedure can 
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be easily implemented whenever required 40]. Regarding this issue we are convinced that, 



despite the present-day lack of consensus among nuclear theorists on how to derive the 
NME's in a direct and controlled manner, they will be able to surmount this obstacle in 
the near future without having to resort to extremely complicated theories. In fact, nuclear 
physics is not merely complicated mathematics: it requires much art to discover the most 
important degrees of freedom and to disentangle the underlying symmetries, for the purpose 
of building very imaginative models, and to skilfully manipulate the pertinent adjustable 
parameters. One should always keep in mind Milton's witty remark: The very essence of 
truth is plainness and brightness. Very likely, once fixed the NME's, it will be possible to 
answer some fundamental questions about neutrinos. 

I'm very sure Dubravko would agree with all that has been said above, and this greatly 
encourages me to pursue the presented line of research and it is a further reason to dedicate 
this article to him. 
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